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ABSTRACT 

Obscured AGNs provide an opportunity to study the material surrounding the central en¬ 
gine. Geometric and physical constraints on the absorber can be deduced from the reprocessed 
AGN emission. In particular, the obscuring gas may reprocess the nuclear X-ray emission pro¬ 
ducing a narrow Fe Ka line and a Compton reflection hump. In recent years, models of the 
X-ray reflection from an obscuring torus have been computed; however, although the reflect¬ 
ing gas may be dusty, the models do not yet take into account the effects of dust on the pre¬ 
dicted spectrum. We study this problem by analyzing two sets of models, with and without the 
presence of dust, using the one dimensional photo-ionization code Cloudy. The calculations 
are performed for a range of column densities (22 < log[A^iy (cm“^)] < 24.5 ) and hydrogen 
densities ( 6 < log[niy (cm“^)] < 8). The calculations show the presence of dust can enhance 
the Fe Ka equivalent width (EW) in the reflected spectrum by factors up to w 8 for Compton 
thick (CT) gas and a typical ISM grain size distribution. The enhancement in EW with respect 
to the reflection continuum is due to the reduction in the reflected continuum intensity caused 
by the anisotropic scattering behaviour of dust grains. This effect will be most relevant for 
reflection from distant, predominately neutral gas, and is a possible explanation for AGNs 
which show a strong Ee Ka EW and a relatively weak reflection continuum. Our results show 
it is an important to take into account dust while modeling the X-ray reflection spectrum, and 
that inferring a CT column density from an observed Ee Ka EW may not always be valid. 
Multi-dimensional models are needed to fully explore the magnitude of the effect. 
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1 INTRODUCTION 


All active galactic nuclei (AGN) are powe red by gas a ccreting 
onto a central supermassive black hole (e.g. fBalbusll2003h . a pro¬ 
cess which emits a significant amount of energy across the elec¬ 
tromagnetic spectrum. Interestingly, a sig nificant number o f AGNs 
show the presence of local obscuration iComastr i||2004|) at dis¬ 
tances ^ 1-10 pc from the central engine (e.g. Antonuccilll993l : 
lUrrv & Padovanil[l995h . The nature and origin of the absorbing gas 
is largely unknown but the gas presents an unique opportunity to 
study material as it transitions from the galaxy to the AGN envi¬ 
ronment. 

Superimp osed onto the t ypical F ^ 1.8-2 X-ray power-law 
in AG Ns fe.g. lDadinall200^: [Beckmann et al.ll2009l: ICorral et al 


20n) is an Fe Kq line (e.g. Nandra & PoundsI 19941 : Ebisawa et al. 


19961) . The narrow c omponent of the lin e is observed in a large 
number of objects l lYaaoob e t al.| |200^: |Kasp il l200lh , even in 
many high-redshift samples terusa et aljboOSi : ICorral et al.ll200^ : 
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IChaudharv et alj |2012| : [l wasawa et alj l2012h . and the width of 
the l ine indicates that it may originate in the obscur ing material 
(e.g. lYaaoob & Padmanabhanl l2004l : IShu et al.l l2010h . Moreover, 
the equivalent width (EW) of the Fe Ka line can be as large as sev¬ 
eral keV fe.g. lLevenson eta lj|2006l) which may be a signature of 
Compton thick fCT) gas (e.g. lMurphv & Yaaoodl2009h . A correla¬ 
tion between the Fe Ka EW and the line-of-sight column density is 
observed for columns > 10^® cm'^ indicating that the line is often 
produced by CT gas in the A GN environment fe.g. lGuainazzi et ^ 
l2005l : lFukazawa et alj2010n . Thus, the Fe Ka line can be used as an 
important proxy to study the properties of the obscuring gas around 
AGNs. 


In recent years, there have been several Monte-Carlo mod¬ 
els of X-ray reprocessing from a CT torus (e.g . Ikedaet ai] l2009l: 
iMurphv & Yagoobll200w [Brightman & Nandra 12011 ah . However, 
all of these calculations omit dust from the models which 
in many cases may be present in th e X-ray absorbing gas 
jjaffe et alJliooi: IPrieto 2004 et alJl20()i: IPrieto 2005 et all 20051: 


iMeisenheimer et al.ll2007l : Tristram et alj|200'^ Raban et al.l2009h . 

As these models are now being used to determine the gross proper- 
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ties of the obscuring material, it is important to consider the effects 
of dust on the X-ray reflection spectrum. Here, we use Cloudy mod¬ 
els to show that the inclusion of dust in the distant obscuring gas 
may have an important impact on the reflected Fe Ka EW. Sect. 2 
describes the geometrical set up of the models and the calculation 
procedure. Sect. 3 describes the results, and Sect. 4 discusses the 
implications of this study. 


2 CALCULATIONS 


The photo-ionization code Cloudy vl3.02 Jperland et al.ll2013b . is 
utilized to study the impact of dust on the reflected iron Ka line. 
Although the calculations with Cloudy are one dimensional, the 
code incl udes multiple grain species with a realistic size distribu¬ 
tion (e.e. lMathis et al.ll l977 ) as well as a compreh e nsive treatment 
of grain physics l Weingartner and Draind 1200 ll : Ivan Hoof et al.l 
|2004 : IWeingartner et al.l l2006h . scattering and absorption of 
dust based on its size and shape. All calculations assume an AGN 
spectrum with an intrinsic luminosity of 10"*^ er g s'* in the 2-10 


keV b and and a typical photon index of 1.9 fe.g. iBeckmann et aH 
1 2009 ). The optical to UV spectrum is characterized by aox = -1.4 
(e.g. Zamoratii et aLlllOsH) . a UV slope of -0.5 fe.g. iFran^l 19931 : 
lElvis et alJll994h . and a blackbody temperature of 1.4 x 10® K (ap¬ 
propriate for a black hole of lO’ Mq accreting at 10% of its Ed¬ 
dington rate). The gas irradiated by this spectrum has an uniform 
hydrogen density, uh, extended radially outward until it reaches a 
specified column density, Nh- The covering factor is fixed at 0.67, 
equivalent to a 30° opening angle, and calculations are performed 
for 22 < log[Vir(cm“^)] < 24.5 (with 0.25 dex spacing), and 
6 < log[nir(cm“®)] < 8 (with a spacing of 1 dex). 

To study the effects of dust, two sets of m odels are com¬ 
puted . The first has gas with solar abundances dOrevesse et al.l 
120101) and n o dust (ND) , whil e the second include s ISM 
dust (WD) I Mathis et alj 19771 : Cowie and Songail^ 19861 : 
Savage and Sambach et alj 1 19961 : Meyer et ^ 19981 : Snow et all 


20071). In order to isolate the effects of dust on the reflected spectra. 


the abundances in the gas in the WD models are kept at solar, but 
the metal abundances are redu ced using the depletion factors of 
ICowie and SongaiT3 lll986l) . and |jenkini ( ll987t) and the gas-to-dust 
ratio was kept at 151. To check the self consistency of the two 
models, the abundances of each element are added and the overall 
abundance of elements is consistent within 0.05% in both ND and 
WD models; however, the Fe abundance in the WD models is 5% 
higher than ND models. 

At high AGN luminosities grains may reach the sublimation 
temperature and return the metals back to the gas phase. This effect 
is not easily treated by Cloudy, so the inner radius of the illuminated 
gas is fixed at a distance of 10 pc away from the AGN. At this 
distance the dust temperature inside the illuminated gas is always 
below the sublimation temperature of the grains. As a further check 
on the choice of parameters, the reprocessed 12/rm luminosity of 
the WD models is compared again st the one expected from the 12 
^m-X-ray luminosity relationship jGandhi et al.l200^ . Depending 
on Uh, the simulated 12/rm luminosities are 1.2-1.5 times larger 
than the observationally predicted one. Thus, our simple dusty gas 
setup will be a reasonable model for the average properties of a 
Seyfert-type AGN. 

Finally, to study the effects of dust on the Fe Ka emission 
line, the equivalent width (EW) of the line is analyzed as a function 
of Nh and uh for both sets of models. The EW is defined as the 
ratio of the intensity of the Fe Ka emission line (IfuKc) to the 


reflected continuum at the line energy (Jc) multiplied by the energy 
bin width (AE): 


EW 


IpeKa 

Ic 


X AE. 


( 1 ) 


Ideally, E should be measured at 6.4 keV, but because of the pres¬ 
ence of the emission line, it is difficult to determine E at that 
energy. Therefore, E is measured at 6.3 keV as there are no pre¬ 
dicted emission lines at that energy. Using a different E between 
6.3 keV and 6.5 keV leads to only very small changes in the EW. 
Rather than measuring IpeKa from the predicted spectrum, we use 
the emergent line intensities reported in the Cloudy output since 
it takes into account the effects of extinction on the line intensi¬ 
ties. In the ND models, the Ee Ka emission line comes from cold 


sas (e.e Georse & Fabianll991 ; Weaver. Gelbord. & Yaaoobl200ll; 

Paee et alj 120041; lYac 

oob & PadmanabhanI |2004|; IZhou & Wand 

20051; Jiane et alj2006 

; Levenson et al.l2006l). and in the WD mod- 


els it comes from both cold gas and grains which are added to ob¬ 
tain IpeKa- The cold Fe Ka fluorescence line defined by Cloudy 
is actually em itted between 6.4 to 6.424 keV for Fe I to Fe XIV 
jHousell 19691) : therefore, we set AE = 24 eV in eq. 1. As a con¬ 
sistency check, the EW calculation was compared with the result 
from Ikeda et al. (2009). Our results predict the ND Fe Ka EW to 
be « 1.8 keV for Nh = 10^“* cm~^ and uh = 10^ cm“® which 
is slightly lower tha n the EW ss 2 keV predicted from the Monte 
Carlo simulations of llkeda et al.l ( l2009l) . However, given the differ¬ 
ences in computational techniques and physical setup, this differ¬ 
ence is adequate. In addition, we are focused on the relative changes 
in the EW when dust is included in the irradiated gas. 


3 RESULTS 

Figure 1 shows how the Fe Ka EW depends on Nh and uh ■ The 
left panel of figure 1, shows that the Fe Ka EW increases with Nh 
for Compton thin gas (10^^ < Nh < lO^"* cm “). Such a depen¬ 
dency is expected since the amount of illuminated gas increases 
with Nh and this produces more Fe Ka emission. For Compton 
thick gas {Nh > 10^"^ cm“^), the amount of gas irradiated by X- 
rays is limited by Compton scattering and the Fe Ka EW is fairly 
invariant with Nh ■ In addition, the EW of Ee Ka increases with 
Uh in both the ND and WD models due to the rise in total opacity 
with Uh- This reduces E and therefore increases the EW. 

Our calculation shows that the presence of dust may signifi¬ 
cantly enhance the Ee Ka EW. The right panel of figure 1 shows 
that the EW of the Ee Ka line is increased by a factor of 2.5 for 
Uh = 10® cm ® (for Nh > 10®'* cm ®) due to presence of grains. 
For Uh = 10® and 10* cm ®, the EW is enhanced by a factor of 
4.1 and 8, respectively. To illustrate the reason for this enhance¬ 
ment, Fig. 2 plots E and IpsKa ratios vs. Nh - This figure clearly 
shows that the change in EW between WD and ND is due to the 
reduction in the continuum and not because of a change in IpeKa- 
{IpeKa does slightly increase due to the presence of grains, but 
this is largely due to the abundances in the models ND and WD not 
being perfectly self consistent; the Fe abundance is 5% higher in 
WD models. Therefore, it is expected to have slightly more Fe Ka 
emission line in the case of WD.) 

Fig. 2 shows the increase in EW of Fe Ka when reflected by 
dusty gas is due to a reduced continuum. This suppressed reflection 
continuum is due to the reduction in backscattering opacity (Ks) in 
dusty gas. Fig. 3 shows how Ks varies inside the column of an il¬ 
luminated gas as a function of depth for WD and ND models for 
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Figure 1. Left : Solid and dashed lines plot the EW versus for the case of ND and WD respectively. The EW increases with N[f and njf in both ND 
and WD cases. Right : The Ee Ka EW gets significantly enhanced in the WD models with a factor that depends strongly on . 



loglNn (cm^)] 



Figure 2. Left : The effect of dust on R- As nn increases, the effects of dust grows and the R ratio significantly decreases and is therefore responsible for 
the increased Fe Ka EW. Right : The presence of grain only slightly increases the Fe Ka intensity, but this has a very small effect on its EW. 


Nh = cm^ and uh = lO’ cm‘^. At depths > 10*^ cm, H is no 
longer ionized and dust contributes significantly to the scattering 
opacity; therefore, Ks decreases significantly deep inside the cloud 
in WD models compared to the ND models. This drop in k® due to 
dust is because scattering from grains is highly anisotropic and fa¬ 
vor scattering in the forward direction by an amount dep ending on 
the size, structure and shape of the grains dPrain 3 l2003li . Since the 
wavelength of X-rays is smaller than the size of most grains, there 
is only a weak coupling between the radiation and the grain lead¬ 
ing to anisotropic scattering. As the Fe EW is computed using the 
backscattered continuum, this anisotropic scattering will reduce the 
backscattered intensity for WD models and thus the EW increases. 


It has been proposed that dust in an AGN enviro nment has 
fewer small grains than the typical ISM di stribution jMaiolind 
l200lh . As ka is proportional to grain size jHavakawal 1970h . a 
grain size distribution deficient in small grains will likely produce 
a stronger reflection continuum and thus the Fe Ka EW will not 
increase as strongly as with the ISM grains. To check this effect we 
computed a Cloudy model (n_tf=10’ cm‘^, Nh = 10^“* cm"^) with 


Orion abundances that includes a grain size distribution deficient 
in small grains. The ka from this model is overplotted on Fig. 3 to 
compare against the results from the ND and WD calculations. As 
expected, the Orion ka is larger than the one using ISM grains, but 
still smaller than the model with no dust. The EW of the Fe Ka line 
is enhanced by a factor of « 4 with the ISM grains, and by a factor 
of « 2.6 with the Orion grains. Hence, the anisotropic scattering 
behaviour of grains may still have an observable impact on the Fe 
Ka EW if the smaller grains are not present. 


4 DISCUSSION AND SUMMARY 

Many groups have computed two dimensional X-ray reflec¬ 
tion models for AGN t orii withou t takin g into account the 


I 2 OO 9 I: [Briehtman & Nandral 1201 la 
sible that distant X-ray refiectoi 

). However, it is 
s arise from dust' 

pos- 

^ gas 

Jjaffe et alj|2004 IPrieto 2004 et alj|2004 IPrieto 2005 et al.ll2005l; 
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Figure 3. The isotropic backscattering opacity in the irradiated gas as a function of depth for the tlh = 10^ cm'^ and cm'^ models. Deep inside 

the column of gas at depths > 10^^ cm, H recombines and the scattering is dominated by the grains which is highly anisotropic and results in lowering the 
backscattering opacity in the WD models. The Orion WD model shows a higher ks than the ISM WD model due to the deficiency of small grains. However, 
the Orion disti'ibution still produces a lower backscattering opacity than the ND model which corresponds to an increase in the Fe Ka EW by a factor of 
2 . 6 . 


iMeisenheimer et alJ20Q7l : iTristram et alJl2007l : iRaban et al.ll2Q09l) . 

Moreover, the Fe Ka line is an impor tant proxy to estimate the 
colum n density of distant observing gas (iGhisellini. Haardt & Matfl 
Il994h . Therefore, we used Cloudy to compute simple one¬ 
dimensional models of dusty gas illuminated by an AGN in order 
to study the effects of dust on Fe Ka emission line. We found that 
the presence of dust may significantly enhance the Fe Ka EW (by 
factors of ~ 5) in the reflection spectrum even in non-CT gas. 
When grains are present in the gas, scattering is anisotropic and 
there are less backscattered photons in the reflected continuum and 
the overall continuum intensity is decreased. This suggests that 
inferring a CT Nh from the Fe Ka EW can be precarious. 

The increase in EW occurs when dust dominates ks and there¬ 
fore will be most important when the gas contains predominantly 
neutral hydrogen. This limits the reflecting cloud to be relatively 
distant from the nucleus or to have a significant density. There¬ 
fore, the Ee Ka EW enhacement may only be important for a 
certain subset of AGNs that exhibit infrared emission from AGN 
heated dust, a large Ee Ka EW, and an unusually weak Comp¬ 
ton reflection component. Eor exa mple, NGC 7213 is observe d to 
have significant hot dust emission dRuschel-Dutra et^l2014l) . an 
Ee K a EW ~ 120 eV, and a very weak Compton reflecti on compo¬ 
nent dSianchi et al.l200^ : lEmmanoulopou los et al. l2013l). Simi larv. 
NGC 2210 also shows strong dust emission dHbnig et al.l201^ and 
is observed t o have an Ee Kg EW o f « 35 — 200 eV from Comp¬ 
ton thin gas dMarinucci et al.ll2015h . Alth ough a broad-line re gion 
origin for the Ee Ka line is a possibility dSianchi et al.ll2003h . our 
results also suggest that reprocessing from the dusty Compton thin 
gas in the absorber may also contribute to the Ee Ka line. 

We conclude that the anisotropic scattering behavior of grains 
is an important mechanism to take into account when modeling 
X-ray reflection from torii since it may have a significant effect 
on the predicted Ee Ka EW. However, the effects of geometry are 
important. Since the scattering behavior of grains directly affects 
the continuum, the effects of grains depends on the viewing angle 
relative to the orientation of torus, therefore the true magnitude of 
the effect needs to be studied in multi-dimensional models. The EW 


of the emission line may be reduced, enhanced or unchanged since 
the continuum can be the transmitted spectrum, reflected spectrum, 
incident spectrum or their combination depending on the viewing 
angle. 
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